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Abstract. Radio galaxy PKS 1138—262 is a massive galaxy at z = 2.16, located in a dense environment. We have 
found an overdensity of Lya emitting galaxies in this field, consistent with a proto-cluster structure associated 
with the radio galaxy. Recently, we have discovered forty candidate Ha emitters by their excess near infrared 
narrow band flux. Here, we present infrared spectroscopy of nine of the brightest candidate Ha emitters. All these 
candidates show an emission line at the expected wavelength. The identification of three of these lines with Ha is 
confirmed by accompanying [N n] emission. The spectra of the other candidates are consistent with Ha emission 
at z ~ 2.15, one being a QSO as indicated by the broadness of its emission line. The velocity dispersion of the 
emitters (360 km s _1 ) is significantly smaller than that of the narrow band filter used for their selection (1600 
km s _1 ). We therefore conclude that the emitters are associated with the radio galaxy. The star formation rates 
(SFRs) deduced from the Ha flux are in the range 6-44 Mq yr _1 and the SFR density observed is 5-10 x higher 
than in the HDF-N at z = 2.23. The properties of the narrow emission lines indicate that the emitters are powered 
by star formation and contain very young (< 100 Myr) stellar populations with moderately high metallicities. 
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Cosmology: observations - Cosmology: early Universe 



1. Introduction 



High redshift clusters are prime subjects for the study 
of galaxy formation and cosmology. The powerful radio 
galaxy PKS 1138—262 at Z = 2.156 appears to be the 
brightest galaxy in a high redshift cluster. We have dis- 
covered an ov erdensity of Lya emitters within 1.5 Mpc 
of 1138-262 l|Kurk et al.ll2000t Paper I). The redshifts 
of 14 emitters were spectro scopically confirmed to be in 
the range 2.14 < z < 2.18 l|Pentericci et al.ll2000l Paper 
II). In addition., we h ave carried out near infrared imag- 
ing ftCurk et alJl2004l Paper III). The number of K band 
galaxies and extremely red objects in this field is higher 
than in blank fields. We found 40 objects with excess nar- 
row band flux, consistent with Ha emission at z ~ 2.16. 
The surface density of Ha emitters increases towards the 
radio galaxy and their average K magnitude is lower and 
therefore their inferred stellar mass higher than for the 
Lya emitters. Here, we present infrared spectroscopy of 
nine candidate Ha emitters to confirm their redshift and 
determine the velocity dispersion of the sample. We as- 
sume a flat Universe with ho = 0.65 and ft m = 0.3. 



2. Selection, observations and data reduction 

With the aim of detecting Ha emitting galaxies in the 
proto-cluster associated with PKS 1138—262, imaging in 
K s and in a narrow band filter (A c = 2.07^, FWHM = 
0.026^) was carried out, employing two pointings cover- 
ing a total field of 12.5 arcmin 2 . There are 40 objects with 
rest frame equivalent width (EWo) > 25 A (see Paper 
III). From the list of 29 candidate Ha emitters within 
1 '3 from the radio galaxy, we selected those with Ha flux 
> 3.5 x 10 -17 erg cm -2 s _1 . Furthermore, we selected those 
that were conveniently located for placement in the slit 
for spectroscopic follow-up, which was carried out with 
ISAAC at VLT Antu (UT1) 1 . The short wavelength cam- 
era of ISAAC is equipped with a Rockwell Hawaii 1024 2 
pixel Hg:Cd:Te array which has a projected pixel scale of 
0.147". We used the medium resolution grating in second 
order resulting in a dispersion of 1.23 A. The observations 
were carried out in the nights of March 23 and 25, 2002 
under variable seeing, which was just below 1" for most 
of the time. The l"xl20" slit employed resulted in a res- 
olution of ~ 2600. During acquisition, the slit was first 
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positioned on a bright point source (a star or the radio 
galaxy) and subsequently positioned at the midpoint be- 
tween two candidates, which was always within 32". We 
have employed four slit positions, each targeted at two or 
three candidate emitters for 3.5 hours (3.75 in one case). 
In total, nine candidates were observed, one of which was 
included in two slits. Per slit, we obtained 14 (15) frames 
of 15 minutes with offsets of 15, 18 or 20" in ABBA se- 
quence with additional random jitter offsets up to 5, 2, or 
1", respectively, to avoid recurrent registration of spectra 
on bad pixels. Standard stars were observed with the same 
slit at a range of airmasses during the nights to correct for 
telluric absorption and to calibrate the data in flux. All 
observations were carried out at airmass below 1.8. 

Standard data reduction was carried out using pair- 
wise frame subtraction, resulting in a final frame with 
positive spectra accompanied by negative spectra on both 
sides. Care was taken during setup to ensure that these 
negative spectra did not overlap with positive ones. One 
dimensional spectra of the candidate emitters were ex- 
tracted from the positive two-dimensional spectra using 
the spatial profile of a standard star spectrum observed 
during that night and averaged per 3 pixels yielding bins 
of 3.7 A. The spectral resolution is 7 A, as measured from 
the FWHM of the skylines. The wavelength calibration 
is based on the OH skylines observed. The telluric stan- 
dards observed at a range of airmasses show only small 
variations (less than a percent on average). An average 
absorption spectrum per night was used to correct for 
telluric absorption. There are no spectrophotomctric flux 
standards in the infrared. We have therefore used one of 
the telluric standards (Hip043868) with spectral type Bl. 
Such stars have a featureless spectrum in this wavelength 
region given by a black body curve at T = 25,500 K. The 
curve was normalized to the K magnitude of the star and 
subsequently used to divide the extracted and absorption 
corrected spectrum to obtain the flux calibration. The two 
dimensional spectra of all objects are shown in Fig. ^ 



3. Results 

All candidates observed show an emission line, which 
means that our selection was 100% efficient. We have fit 
Gaussian curves to the emission lines applying a least 
squares method (see Fig. |2J in order to determine their 
central wavelength, deconvolved FWHM (if resolved) , flux 
and EWo (Table ^) . Also presented in this table is the 
SFR derived from the Ha and UV luminosities (see dis- 
cussion in Paper III), which is in the range 6-44 M Q 
yr . We do not detect continuum emission in most of 
the spectra, except in the co-added spectrum of candidate 
131, 215 and the serendipitous object in slit one (0.6±0.6, 
0.7±0.6 and 1.7±0.7 xl0~ 18 ergcm" 2 s" 1 A -1 , respec- 
tively) , which compare favourably with the line subtracted 
broad band fluxes measured by imaging in paper III (0.6, 
1.0 and 2.1 xlO" 18 erg cm" 2 s" 1 A -1 , respectively). The 
EWq is therefore based on the line flux measured in the 



spectra and the broad band magnitude measured on the 
images. 

The spectra of objects 131, 183 and 229 were fit by two 
Gaussians for which the relative centers were fixed as for 
Ha and [Nn]A6583A. Candidate 329 shows two emission 
peaks only 17.8±2.7 A apart. Both lines could be Ha emis- 
sion from one galaxy with two components, seperated by 
~ 250 kms" 1 in velocity. We have considered the possibil- 
ity that the lines are due to [On] emission at z = 4.5674 
which would have a separation of 15.0 A. The emission 
line ratio of the supposed [On] A3729/3726 would be 0.7 
which implies an electron density ~ 10 3 cm" 3 , normally 
only observed in the central parts of nebulae ijOsterbrockl 
1989). In addition, a faint (B = 27.0) counterpart in the 
B band, sampling a wavelength range below 912 A for 
z = 4.6, makes the identification with [On] improbable. 
Candidate 131 was included in two slits, with a position 
angle difference of 40°. The two fits to the spectra of emit- 
ter 131 indicate a velocity difference of 180 kms" 1 which 
can be explained by the fact that different regions of the 
galaxy have been sampled. 

For how many of the emission lines can we be sure 
that the identification with Ha at z ~ 2.15 is correct? 
For the three objects with confirming [Nil] lines, we can 
be certain. For the QSO object 215, [Nil] is blended with 
the very broad Ha line and impossible to discern. Given 
a [Nil] /Ha line ratio of 1/3, we do not expect to detect 
the [Nil] line for objects 79, 207, 284, 329 above the noise 
and the identification with Ha is therefore consistent but 
not 100% certain. An identification with Ha for object 144 
can only be true if the [Nil] /Ha ratio is < 1/6 (object 183 
has an observed ratio of 1/5.5). An alternative identifica- 
tion with [O in] A5007A is improbable, as we do not detect 
its counterpart [Olil]A4958A at 2.045/x. We consider this 
therefore a probable Ha identification. 

At least three of the emission lines are spatially ex- 
tended. In two cases, we detect unordered velocity struc- 
ture, but the morphology of the two dimensional spectrum 
of candidate 229 resembles a rotation curve. A fit to this 
structure results in a rotation velocity of 50 kms" 1 at 6 
kpc radius, implying a dynamical mass of 3.5 x 10 9 Mg. 

4. Discussion 

In Fig. |3 we show the redshift distribution of the 9 emit- 
ters, assuming that all emission lines can be identified as 
Ha. For this plot, we used the average redshift of the two 
lines of object 329. Also plotted is the sensitivity curve 
of the narrow band filter used to select the candidates. 
A random distribution of emitter redshifts would follow 
this curve. We have run a Monte Carlo simulation with 
10,000 realizations of a randomly sampled distribution of 
emitters given this filter as selection criterium. Although 
the mean of the measured distribution is only 0.18<r away 
from the mean of a random sample, its dispersion is much 
smaller, deviating by 1.75<r from a random sample. The 
probability that the redshifts we measure are drawn from a 
random distribution is therefore 8%. Note that the redshift 
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Fig. 1. The ten two dimensional spectra observed through four slits. Skyline residuals are visible as vertical lines with 
a higher noise level. Slit and spectrum number are indicated, as are the object number from the Ha candidate list. The 
images have been smoothed with a mask of 3x3 pixels. Near the emission line of candidate 229, continuum emission 
of a serendipitous galaxy is visible. The spectrum of candidate 215 shows a very broad line which almost covers the 
full spectral range. 



Table 1. Properties of the observed emission lines 



ID 


S 


R.A. & Dec 


. (,12000) 


Line 


Fit 




2 


FWHM 


Flux 


EW 




Ha SFRuV 


O/H 


(1) 


(2) 


(3) 




(4) 


(5,6) 




(7) 


(8) 


(9) 


(10) 




(11) 


(12) 


(13) 


79 


4 


11:40:52.62 - 


26:30:01.0 


Ha 


1,2 


2 


1558 


160± 80 


3.6±1.9 


90± 


50 


11±6 


33±1 




131 


2 


11:40:51.28 - 


26:29:38.7 


Ha 


3.1 


2 


1518 


320± 50 


12.3±2.2 


65± 


15 


39±7 


19±1 


8.8 










[Nil] 








270±100 


4.4±2.0 


25± 


10 






±0.2 


131 


4 


11:40:51.28 - 


26:29:38.7 


Ha 


3,1 


2 


1537 


360± 80 


8.2±2.4 


45± 


15 


26±8 


19±1 


8.8 










[Nil] 








420±250 


3.4±2.5 


20± 


15 






±0.2 


144 


3 


11:40:43.45 - 


26:29:37.5 


Ha 


1,3 


2 


1463 


40± 20 


6.7±1.3 


330± 


70 


21±4 


15±1 




183 


2 


11:40:46.15 - 


26:29:24.9 


Ha 


3,1 


2 


1546 


170± 20 


13.8±1.7 


230± 


35 


44±5 


40±1 


8.6 










[Nil] 








200±120 


2.5±1.7 


40± 


30 






±0.4 


207 


4 


11:40:50.20 - 


26:29:21.0 


Ha 


1.2 


2 


1540 


<100 


1.9±0.9 


50± 


25 


6±3 


12±1 




215 


1 


11:40:46.01 - 


26:29:16.9 


Ha 


3,1 


2 


1568 


5300±800 


46.2±8.8 


150± 


30 


t 






229 


1 


11:40:46.10 - 


26:29:11.5 


Ha 


3,1 


2 


1489 


290± 60 


7.1±1.9 


30± 


10 


23±6 


26±1 


8.8 










[Nil] 








130± 90 


2.4±1.7 


11± 


8 






±0.3 


284 


3 


11:40:45.58 - 


26:29:02.4 


Ha 


1,2 


2 


1523 


90± 50 


3.0±1.4 


300±140 


10±4 


6±1 




329 


3 


11:40:46.88 - 


26:28:41.4 


Ha 


2.3 


2 


1609 


60± 40 


2.7±1.3 


1350±700 


9±4 


7±1 












Ha 


2,3 


2 


1636 


110± 30 


2.0±1.1 


1000±550 


6±4 







Notes: (1) Candidate number (2) Slit number (3) Coordinates (4) Line identification (5) Type of fit: 1) one Gaussian curve, 
2) two Gaussian curves, 3) two Gaussians for Ha and [Nil] (6) Identification is 1) certain, detection of [Nil], 2) consistent, 
expected non detection of [Nil], 3) possible (7) Redshift with random error 0.0002, except for candidate 215, for which it is 
0.002 (8) FWHM in km s" 1 (9) Flux in 10~ 17 erg cm" 2 s" 1 (10) EW in A (11) SFR Ha in M Q yr" 1 (12) SFRuv in M yr _1 
(13) MetaUicity in 12 ± log[0/H] ( f ) Ha unrelated to SFR. 

of the radio galaxy is 2.156 l|Rottgering et al.lll997th very 
close to the mean of the selection filter (2.152), and to the 
mean of the measured redshift distribution (2.154). The 
distribution is consistent with a group of Hcv emitters as- 
sociated with the radio galaxy. The velocity dispersion of 
this group is 360 kms -1 (using the gapper sigma method, 



Ifieers et alFfl OOL while the virial radius of the nine emit- 
ters is 0.45 Mpc, implying a virial mass of 8 x 10 13 M Q (as- 
suming that all lines can be identified with Ha and taking 
the mean of the two redshifts for object 131). This mass 
is merely illustrative as at this redshift it is improbable 
that the structure is virialized. The velocity dispersion of 
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2.060 2.065 2.070 2.075 2.060 2.065 2.070 2.075 2.060 2.065 2.070 2.075 

Fig. 2. One dimensional spectra of the eight narrow line 
emitters (histograms) with fits overlayed (solid lines). 
Object 329 has also an alternative fit for [Oil]. Units are 
in 10 -18 ergcm~ 2 s _1 A -1 and fim. The grey band denotes 
the position of a sky line. The dashed line indicates the 
position of the [N II] line for identification of the detected 
line with Ha. Note the deviant flux density scale in the 
last row. 



the Ha emitters is smaller than the velocity dispersion of 
the confirmed Lya emitters, both for the complete sample 
(1050 kms -1 ) and the nine within the solid angle of the 
two ISAAC fields (760 kms" 1 , shaded part of histogram 
in Fig. |2l ■ There is no evidence for a bimodal redshift 
distribution as observed for the Lya emitters. 

We can construct a complete sample out of the spec- 
troscopic sample by excluding the two objects with the 
lowest Ha flux and including the radio galaxy. This col- 
lection represents all candidate Ha emitters with Fjja > 
4.0 x 10~ 17 ergcm~ 2 s -1 within 1 '3 from the radio galaxy. 
The FWHM of the narrow band filter (2.134 < z < 2.174) 
and the solid angle given above define a comoving volume 
of 815 Mpc 3 , resulting in a volume density of 0.010 Mpc -3 , 
which is a factor four higher than the density of confirmed 
Lya emitters in this field. All star forming objects detected 
have line fluxes lower than the high redshift Ha surveys 
discussed in Paper III, but we can compare the SFR den- 
sity to the density at z — 2.23 d erived from Ha emis - 
sion in the HDF-N as measured bv llwamuro et al. (2000). 
Following their cosmology (ho = 0.5, qo = 0.5) and pro- 
cedure to correct for the part of the Ha luminosity func- 
tion below the detection limit, we obtain a SFR density 
of 0.48 M yr' 1 Mpc~ 3 . This is 10 (5) times higher than 
the (reddening corrected) value obtained for the HDF- 
N. Using the redshift range defined by the Ha emitters 
(2.146 < z < 2.164) results in values that are larger by a 




2.12 



2.14 



2.16 



2.18 



2.20 



Fig. 3. Redshift histogram for the Ha emitters (solid) and 
the previously known Lya emitters (dotted, and shaded 
for those within ISAAC field). An arrow denotes the red- 
shift of 1138—262. Also shown is the sensitivity curve of 
the narrow band filter used for the selection of the Ha 
candidates. The surface beneath the curve is normalized 
to nine. 

factor of two. Likewise, smaller SFR densities would result 
if some of the Ha lines have been misidentified. 

The properties of the detected emission lines pro- 
vide information about the physical conditions in the 
galaxies. The FWHM of the narrow nebular emission 
lines detected are in the range 40 - 360 kms -1 with 
an average of 190 kms -1 . These values are compara- 
ble to those found for LBGs at z ~ 3 by IPettini et alJ 
l|200ll) . The [Nil] /Ha ratio can be used to distinguish 
narrow-line active galaxies f rom Hn region- like galaxies 
l|Veilleux fc Osterbrockll98"7j) . The three emitters with de- 
tected [Nil] have -0.74 < log([Nn]/Ha) < -0.42, which 
puts them among the star forming galaxies. In the ab- 
sence of shock excitation, the [Nil] /Ha ratio can also be 
used as metal licity indicator. Using the empirical relation 
calibrated bv lDenicolo et afl ((2002) , the average ratio of 
the three emitters implies 12 + log(0/H) ks 8.7. This 
value is comparable to the broa d range of values obt ained 
for present-day spiral galaxies l(van Zee et al.lfl998h . The 
EWo of some detected narrow lines are surprisingly high, 
up to 1350 A. This can be explained by very young stellar 
populations where the continuum radiation around 6000 A 
is still very weak. EWo values between 200 an d 330 A im- 
ply an age < 100 Myr l|Leitherer et alJll999h . The mod- 
erately high metallicities found for the objects with de- 
tected [Nil] emission, however, require that the galaxies 
are near the end of the star formation event. This require- 
ment seems to indicate that these emitters have undergone 
a very similar evolution. 

5. Conclusion 

Infrared spectroscopy has established the presence of nine 
line emitters within 0.6 Mpc of the HzRG PKS 1138-262. 



J. D. Kurk et al.: A search for clusters at high redshift 



5 



Three emitters show an additional line which confirms the 
identification with Ha at z = 2.15, while four more have 
spectra consistent with Ha at this redshift, one being a 
QSO as indicated by the broadness of its emission line. 
One emitter shows only a single strong line, which is pos- 
sibly Ha and one emitter exhibits two lines which probably 
originate from two emission line regions within one galaxy 
at z — 2.16. Additional evidence for identification of all 
observed lines with Ha is the small velocity dispersion 
(360 km s _1 ) as compared with the width of the selection 
filter. This dispersion is also smaller than the dispersion 
of the Lya emitters. The star formation rate density of 
the observed emitters is a factor ten higher than found 
at z = 2.23 in the HDF-N. These results support the for- 
merly advocated ideas that PKS 1138—262 is located in 
a proto-cluster at z = 2.16. The properties of the narrow 
emission lines indicate that the emitters are powered by 
star formation and contain very young (< 100 Myr) stellar 
populations with moderately high metallicities. It seems 
that we observe these galaxies near the end of their first 
and major burst of star formation. 
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